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(Nj ■ ABSTRACT 



Multifrequency polarimetry with the VLBA confirms the previously reported time- 
varying Faraday rotation measure (RM) in the quasar 3C 279. Variabihty in the RM 
and electric vector position angle (EVPA) of the jet component (C4) is seen making it 
an unreliable absolute EVPA calibrator. 3C 273 is also shown to vary its RM structure 
^ ■ on 1.5 year time-scales. Variation in the RM properties of quasars may result from a 

Faraday screen which changes on time-scales of a few years, or from the motion of jet 
components which sample spatial variations in the screen. A new component emerging 
from the core of 3C 279 appears to be starting to sample such a spatial variation. Future 
monitoring of this component and its RM properties is suggested as a diagnostic of the 
narrow line region in 3C 279. We also present a new method of EVPA calibration using 
' the VLA Monitoring Program. 
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1. Introduction 

The parsec-scale cores and jets of quasars have typical linear polarizations of 1-10% (Cawthorne 
et al. 1993) with the jets being more strongly polarized. If a change in the electric vector position 
angle (EVPA) of the polarized emission is seen with frequency, and can be attributed to the Faraday 
effect, then we can learn about the densities and magnetic fields of thermal gas along the line of 
sight. 

High angular resolution polarimetry has shown that some quasars have Faraday Rotation 
Measures (RM) in excess of ~ 1000 rad m~^. In the case of the quasar 3C 279 these high RMs 
can change over 1.5 years by over 1000 rad m^^ (Taylor (2000) and references therein). These 
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high RMs are present within a projected distance of ~ 20 pc from the core of the quasars. Beyond 
this distance the RMs decrease to < 100 rad m~^, a value consistent with that expected from the 
passage of the emission through the host galaxy and through the ISM of our own Galaxy. 

In §2 we describe VLBA observations of the quasars 3C 273 and 3C 279, with the results of 
the rotation measure analysis presented in §3. In §4 we discuss our results in light of the unified 
scheme of AGNs. 

We assume Hq = 50 km s~^ Mpc~^ and qo=0.5 throughout. 

2. Observations and Data Reduction 

The observations, performed on 2000 January 27 (2000.07), were carried out at ten frequencies 
between 8.1 and 43.2 GHz using the 10 clement VLBA^. Sources 3C273, 3C 279, and calibrator 
J1310+3220 were observed for approximately 2.5 hours each. The calibrator OJ287 was observed 
for almost one hour. Right- and left-circular polarizations were recorded using 1 bit sampling 
across a bandwidth of 8 MHz. The VLBA correlator produced 16 frequency channels across each 
IF during every 2 s integration. 

Amplitude calibration and fringe-fitting were performed as in Taylor (2000). Weather, in 
particular snow on the antenna, caused us to flag the data from Hancock at 43 GHz. We also 

encountered radio-frequency interference at 12 GHz on our 3C 279 scans, which we attribute to 
geostationary television broadcasts. This required extensive editing of the 12 GHz data on 3C 279. 

Feed polarizations of the antennas were determined using 3C 279 and the AIPS task LPCAL. 
We assumed that the VLBA antennas had good-quality feeds with relatively pure polarizations, 
which allowed us to use a linearized model to fit the feed polarizations. Absolute electric vector 
position angle (EVPA) calibration at 8-22 GHz was determined by using the EVPA of OJ287 
obtained on 2000 January 29 listed in the VLA Monitoring Program^ (Taylor k. Myers 2000). 
Polarization monitoring observations of J1310-t-3220 on 1999 December 19 and 2000 February 18 
were interpolated to the observation date and employed as a check of the O J287 calibration (Table 
1). At 43 GHz the low polarized flux density from these calibrators prevented an absolute EVPA 
calibration. To set the EVPA at 43 GHz we performed a linear least-squares fit to the 3C 279 
component C4 EVPA data, and used this fit to set the EVPA for 3C 279 at 43 GHz. The correction 
applied for 3C 279 C4 was then used for the other sources (Fig 1). 

Taylor (1998) suggested that C4 may be used as an absolute EVPA calibrator due to its low 
RM and strong polarized and total flux densities. From 1994 (Leppanen et al. 1995) through 1998 
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(Lister et al. 1998; Wardle et al. 1998; Denn et al. 2000; Taylor 2000) component C4 maintained 
an EVPA oriented nearly parallel to the jet axis. If we define the intrinsic polarization angle (IPA) 
of a source as the intercept in the RM plot (Simard-Normandin et al. 1981) the IPA of 3C 279 C4 
determined by our epoch 2000.07 observations is —87.1 it 2.2°. The position angle of component 
C4 is significantly different from the value of ~ —120° ±3° found by Taylor (1998). Therefore 
we conclude that C4 cannot be used as an absolute EVPA calibrator, without complimentary 
monitoring observations. 

3. Results 

We present multifrequcncy radio observations, polarization properties, and rotation measure 
distributions for 3C 273 and 3C 279. In the rotation measure images presented in figures 2 through 
5 the higher frequencies have been tapered to provide a matched restoring beam of the dimensions 
indicated in the captions. 

3.1. 3C273 

Figure 2 presents a time series RM mosaic for this quasar (epoch 1997.07 from Taylor 1998). 
Beyond 10 pc from the nuclear region the jet exhibits a relatively smooth and uniform RM distri- 
bution of ~ +500 rad m~^. The RM of this region remains constant over a 3 year time span. This 
jet continues to show good agreement with the law of Faraday rotation. In contrast the nuclear 
region changes from an apparently smooth —2000 rad m~^ in 1997 to a more complex structure 
which varies by | 4000 | rad m^^ across 10 pc. To investigate the structure of this complex central 
region we created a higher resolution RM map in Figure 3. This map consists of data from 15-43 
GHz with the high frequency data tapered to the resolution of the 15 GHz data. Inserted in this 
figure are the data used to obtain the indicated rotation measures. All points used to calculate the 
RM are at least 1% polarized. The central region consists of sub-components of RM in the range 
I 1000 — 2000 I rad m^^. Good agreement to a law is present except in the northeastern-most 
part of the source, closest to the center of activity. 

3.2. 3C279 

Figure 4 presents a multi-epoch RM map of 3C 279 (first two epochs from Taylor 1998 and 
Taylor 2000 respectively). The core may be ejecting a new jet component, as shown by the increasing 
extension of the emission to the southwest over the three epochs. Taylor (2000) previously showed 
that 3C279 shows variation in its RM properties. The figure illustrates these changes as the 
central region of 3C279 changes from -1300 to > -500 rad m'^ between the 1997.07 and 1998.5 
observations. 3C 279 maintains this lower RM during our 2000.07 observation with the exception of 
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a "bridge" of high RM which appears to connect the central region with component C4 (Leppanen 
et al. 1995). 

The high resolution RM structure of this quasar is presented in Figure 5. The rotation measures 
were determined using the 12-43 GHz data tapered to the 12 GHz resolution. The bridge appearing 
in the low resolution data is resolved and two edges of ~ 500 rad m^^ bordering the central region 
of the quasar and C4 appear. Points used to fit the RM in the center of the image and in component 
C4 are at least 1% polarized. Points used in the bridge range from 0.3-0.5% polarized intensity 
and should be interpreted with caution. The asymmetry of the central region is consistent with the 
emergence of a new component from the core. This high resolution image also shows the slope of 
the RM changing from negative at the central peak of emission to a positive slope at the proposed 
newly emerging component. 



4. Discussion 

Taylor (1998, 2000) showed variation in the RM structure of AGNs on small spatial scales. For 
3C279 Taylor (2000) presented time variation of the RM structure of 3C279 over 1.5 years. Our 
2000.07 observations show temporal variations in 3C 273 and demonstrate that this phenomenon 
is not unique to 3C 279. We find good agreement to a law except near the core of 3C 273 
and the resolved bridge of 3C279. The high RMs observed in the central regions could result 
from the passage of the radiation through the narrow-line region (NLR) of the quasar. If the 
measured RM results from a polarized component which dominates within the telescope beam then 
the RM of the central region should vary as new components are created and ejected. As these 
new components emerge from the core their radiation passes through the NLR and a high RM 
results. The component progresses along the jet and the RM varies as the line of sight traverses 
an inhomogeneous NLR Faraday screen. Eventually RMs below 100 rad m^^ should be observed, 
which are comparable to RMs produced by traversing the ISM of the galaxy, as the NLR is exited. 
With this picture the jet of 3C 273 may still probe the NLR as the RM in the jet of this quasar 
in the 8-15 GHz mosaic is approximately 5 times that observed in radio galaxy hosts (Feretti et 
al. 1999). We intend to monitor the cores of quasars for time variations in the RM created as new 
components are ejected and sample these regions. 

A crude estimate of the magnetic field strengths in these quasars can be made by using the 
variance in the observed RMs (Felten 1996). This attributes the observed RMs to a randomly 
ordered magnetic field embedded in cells of uniform size and field strength. From our high resolution 
RM images we estimate a cell size (/) of 1 pc. We estimate a in these maps as 1160 radm^^ for 
3C 273 and 255 radm^^ for 3C 279. With an electron density rig = 10"^ cm^'^, a volume filling factor 
e = 0.01 and a path length X through the NLR of 100 pc the B field can be estimated using 

B = (3/0^/^(a/812eneX). (1) 



This results in a S field of 0.3 mG in 3C 273 and 0.05 mG for 3C 279. This is less than that 



-5- 



estimated by assuming the magnetic and thermal pressures to be in equilibrium for the same 
physical conditions (Udomprasert et al. 1997). This estimate is also less than that estimated for 
the accretion disks of other AGNs obtained by equating magnetic to thermal pressures (Jones et al. 
2000). Direct B field measurements for other galaxies may be useful in setting a scale of reasonable 
B field values to use in interpreting these quasar RM results. Zeeman splitting measurements in 
the Seyfert galaxy NGC 4258 provide an upper limit on the B field within 0.2 pc of the central 
engine of this AGN of 300 mG (Herrnstein ct al. 1998). In the more benign region of the center 
of our own galaxy Zeeman measurements yield B fields varying from 0.5 mG in Sgr B2 ('^ 120 pc 
from the galactic center) (Crutcher et al. 1996) to 4 mG within 2.5 pc from Sgr A* (Plantc ct al. 
1995). It seems reasonable that the field strengths within the NLR of a powerful quasar should be 
at least as strong as that in the center of our own galaxy, so the variance method of Felten may 
provide a lower limit on the B field present in the NLR of these quasars. 

5. Conclusions 

Both 3C 273 and 3C 279 exhibit RM variations on spatial scales of a few parsecs and 1-3 year 
temporal scales. The low resolution RM structure of 3C273 is resolved into sub-components whose 
general structure is visible at low resolution. Good agreement to a law is found outside of 
the core of this quasar. We estimate a B field of 0.3 mG for 3C 273 and 0.05 niG for 3C 279 if 
the RMs result from a tangled field in 1 pc cells. Component C4 of 3C 279 is shown to have a 
varying absolute EVPA and thus is not suitable as a stand-alone EVPA calibrator. We find that 
the VLA Monitoring Program database provides a method of absolute EVPA calibration for VLBA 
observations (Taylor &; Myers 2000). The RM slope change from negative to positive in the new 
component of 3C 279 emerging from the core may also indicate that a new positive RM feature in 
the NLR is being sampled, and its progress will be examined in future VLBA observations. 

This research has made use of the NASA/IPAC Extragalactic Database (NED) which is oper- 
ated by the Jet Propulsion Laboratory, Caltech, under contract with NASA. This research has made 
use of NASA's Astrophysics Data System Abstract Service. R.T.Z. gratefully acknowledges support 
from the NRAO Summer Graduate Student Research Program and NSF Grant HRD-9628730. 
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TABLE 1 
Polarization Angle Calibration 



Source 


u 


•^VLA 


PVLA 


Xvla 


S'vlba 


A^LBA 


XVLBA 


(1) 


(2) 


(3) 


(4) 


(5) 


((-i) 


(7) 


(8) 


OJ287 


8.5 


2.6 


148 


-52.5 


2.5 


140 


26 




22 


2.5 


100 


-53 


1.9 


81 


-61 




43 


2.1 


89 


-37 


1.1 


8 


-88 


J1310+3220 


8.5 


1.7 


56 


-20 


1.8 


52 


58 




22 


1.2 


59 


-27.7 


1.0 


56 


-26 




43 


1.0 


52 


-29 


1.1 


36 


-14 



Notes to Table 1 

Col.(l). — Source name. Col. (2). — Observing band in GHz. Col. (3). — Integrated VLA flux density 
in Jy. Col. (4). — Integrated VLA polarized flux density in mJy. Col. (5). — VLA polarization angle 
(E- vector) in degrees. Col. (6). — Integrated VLBA flux density in Jy. Col. (7). — Integrated VLBA 
polarized flux density in mJy. Col. (8). — VLBA polarization angle (E- vector) in degrees. 
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Fig. 1.— Plot of EVPA vs. for calibrators J1310+3220 (circles) and OJ287 (triangles) and 
component C4 of 3C279 (squares). Filled symbols are EVPA after applying correction based 
on VLA polarization monitoring data for 8-22GHz. Open symbols represent VLA polarization 
monitoring data for the calibrators. The least squares fit shown (dashed line) for C4 was used to 
set the EVPA at 43 GHz for all soTirccs. Uncertainty in position angle is estimated as ±5° for all 
frequencies after EVPA calibration. Slope of the fit for 3C279-C4 equals a RM of 37 rad m~^. 
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Fig. 2. — Rotation measure mosaic of 3C 273 from data at 8 to 15 GHz for the indicated epochs, 
with contours of total intensity at 8 GHz overlaid. The restoring beam has dimensions 1.0 x 2.5 
milliarcsec at position angle 0°. The colorbar ranges from —4500 to +3000 rad m^^. Contours 
start at 30 mJy/beam and increase by factors of 2. 
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Fig. 3. — Rotation measure image of 3C 273 from 2000.07 data at 43 to 15 GHz, with contours of 
total intensity at 22 GHz overlaid. The insets display rotation measure fits at the points indicated 
in rad m~^. The restoring beam has dimensions of 0.6 x 1.2 milliarcsec at position angle 0°. The 
colorbar ranges from —4400 to +2650 rad m~^. Contours start at 20 mJy/beam and increase by 
factors of 2. 
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Fig. 4. — Rotation measure mosaic of 3C 279 from data at 8 to 15 GHz, for the indicated epochs, 
with contours of total intensity at 15 GHz overlaid. The restoring beam has dimensions 1.0 x 2.5 
milliarcsec at position angle 0°. The colorbar ranges from —1310 to +1000 rad m^^. Contours 
start at 25 mJy/beam and increase by factors of 2. 
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Fig. 5. — Rotation measure image of 3C 279 from data at 12 to 43 GHz, with contours of total 
intensity at 22 GHz tapered to 12 GHz resolution overlaid. The insets display rotation measure fits 
at the points indicated in rad m~^. The restoring beam has dimensions 0.7 x 1.7 milliarcsec with 
bpa 0°. The colorbar ranges from —1000 to +1000 rad m~^. Contours start at 30 mJy/beam and 
increase by factors of 2. 



